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We perform a fully nonlinear analysis of superhorizon perturbation in Hoi^ava-Lifshitz gravity, 
based on the gradient expansion method. We present a concrete expression for the solution of 
gravity equations up to the second order in the gradient expansion, and prove that the solution can 
be extended to any order. The result provides yet another example for analogue of the Vainshtein 
effect: the nonlinear solution is regular in the limit A ^ 1 and recovers general relativity coupled to 
dark matter at low energy. Finally, we propose a definition of nonlinear curvature perturbation 7?. in 
Hofava-Lifshitz gravity and show that it is conserved up to the first order in the gradient expansion. 

o 

(N 

CD 

' The new theory of gravitation proposed recently by Hof ava P, [2 lis expected to be renormalizable and unitary. For 
this reason, it has been attracting significant amount of attention 3|-|6| . The theory is power-counting renormahzable 
because of the anisotropic scaling in the ultraviolet (UV), 

Q^. t — ^ X — ^ bx, (1) 

CD 

i-S^ ' with the dynamical critical exponent z > 3. Since this scaling is called Lifshitz scaling, the theory is often called 
Hofava-Lifshitz gravity. 

^SJ Because of the anisotropic scaling ([1]), the time and the space in this theory must be treated separately. Thus, we 

^ , must abandon the 4-dimensional diffeomorphism invariance. Instead, the fundamental symmetry of the theory is the 

' invariance under the so-called foliation-preserving diffeomorphism: 
\ 

g. t^t'{t), x^x{t,x), (2) 

\ which preserves the way the spacetime is foliated by constant-time hypersurfaces. Since this symmetry is smaller than 

■ the counterpart of general relativity i.e. the 4-dimensional diffeomorphism invariance, the structure of the action is 

' less restrictive and allows more parameters. For example, the kinetic part of the action is a linear combination of 

" ^ [ and K'^^Kij with arbitrary coefficients, where Kij is the extrinsic curvature of the constant-time hypersurface and 

J> K = K^^. Hence, the kinetic part can be written as 



Mil 



j Ndt y^(fx {K'^Kij - XK^) , (3) 



where Mpi is the reduced Planck mass and A is an arbitrary parameter. While in general relativity the value of 
A is fixed to unity due to the 4-dimensional diffeomorphism invariance, in Hoi^ava-Lifshitz gravity any value of A is 
consistent with the foliation-preserving diffeomorphism invariance. 

In Hofava-Lifshitz gravity, in addition to the usual tensor gravitons, there is an extra physical degree of freedom 
called the scalar graviton [ll]. In order for the scalar graviton not to be a ghost, the regime 1/3 < A < 1 should 
be excluded [3]. Outside this forbidden interval, the scalar graviton has negative sound speed squared. Therefore, 
in order for the theory to be observationally viable, we need to impose the condition under which the associated 
long-distance instability does not show up [2|. The condition essentially says that A must be sufficiently close to 1 
in the infrared (IR) , and should be considered as a phenomenological constraint on properties of the renormalization 
group (RG) flow. Since the value of A continuously changes by the RG flow, only the regime A > 1 is allowed. 
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When A is very close to 1, the scalar graviton gets strongly coupled If we adopt the usual metric perturbation 
method then we find that higher order terms in the time kinetic part of the action for the scalar graviton become 
larger. This indicates the breakdown of the perturbative expansion in the scalar graviton sector. Note that this 
does not necessarily imply the loss of predictability since all coefficients of infinite number of terms can be written 
in terms of finite parameters in the action if the theory is renormalizable. However, because of the breakdown of 
the perturbative expansion, we need to employ a more or less non-perturbative method to analyze the fate of the 
scalar graviton in the limit A — )■ 1. Such an analysis was performed in for spherically symmetric, static, vacuum 
configurations and it was shown that the limit is continuous and recovers general relativity. This may be considered 
as an analogue of the Vainshtein effect @ . 

The main purpose of this paper is to provide yet another example of the analogue of the Vainshtein effect in Hof ava- 
Lifshitz gravity. For this purpose, we perform a fully nonlinear analysis of superhorizon cosmological perturbation, 
adopting the so-called gradient expansion method [7|. The result is obviously continuous in the limit A — > 1 and 
recovers general relativity coupled to dark matter. 

We also propose a definition of nonlinear curvature perturbation TZ in Hofava-Lifshitz gravity and show that it is 
conserved up to the first order in the gradient expansion. 

The paper is organized as follows. In §|TT1 we briefiy review the basic equations in Hofava-Lifshitz gravity. In IIIIl 
we introduce the gradient expansion method in this theory and present the solution to equations of motion. In § IIVI 
we proposed a definition of nonlinear curvature perturbation TZ and show that it is conserved up to the first order 
in the gradient expansion. § |V] is devoted to a summary of this paper. In appendix we prove that the solution 
obtained in § IIIIl satisfies the momentum constraint in any order of the gradient expansion. 



II. BASIC EQUATIONS 

In this section we review the basic equations of Hofava-Lifshitz gravity, following the notation in [2] , and reformulate 
them in a way suitable for gradient expansion. Basic quantities of Hofava-Lifshitz gravity are the lapse N{t), the shift 
N^{t,x) and the 3-dimensional spatial metric gij{t,x). Combining these quantities, we can construct 4-dimensional 
spacetime metric in the ADM form as 

ds^ = -N^dt^ + g,j{dx' + N'dt){dx^ + N^dt). (4) 

The fundamental symmetry of the theory is the invariance under the so called foliation preserving diffeomorphism: 

t-)-t'{t), (5) 

which preserves the way the spacetime is foliated by constant-time hypersurfaces. 

By requiring invariance under spatial parity and time reflection, the gravitational action is specified as 



^f7 - 



J Ndt y/gd^x {K^■jK'^ - -2A + R + L^>i) 



(6) 



where y/g is the determinant of gij, Kij is the extrinsic curvature defined as 



(7) 



K (= g'^^Kij) is the trace of Kij, is the spatial covariant derivative compatible with gij, and R is the Ricci 
scalar constructed from gij. To lower and raise an index, gij and its inverse g^^ are used. We impose the so called 
projectability condition, namely we require that the lapse function should depend only on time. In order to realize 
the power-counting renormalizability, the higher curvature Lagrangian Lz>i should include up to sixth or higher 
spatial derivatives. For our analysis in this paper, we do not need to specify the concrete form of Lz>i. Moreover, 
for simplicity we shall not include matter action and analyze the pure gravity described by the action ([B]). 
The equation of motion for gij is Sgij = 0, where 



'^9V — 



9ikgjl 



M 



PI 



2 Sig 

Ny/g Sgu 



1 



-Kp,j + 2Kfpkj + -g^jK^'pki - Ag,j - dj 



(8) 
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Here, pij = Kij — \Kgij, £z>iij is the contribution from Lz>i and Gy is Einstein tensor of gij. The trace part and 
traceless part of this equation are, respectively. 



(3A - 1) ( d^K + ^kA + ^A^^A^^ + Z = 0, (9) 



and 

1 



d^A% + KA\ + —iA^^d^N' - A\d,N^) - Z\ - -Z5) = 0, (10) 



where 

3 



EE XV - (11) 



d^^Udt-N^du). (12) 



is the traceless part of if'^-, 

r), = 

N 
and 

ZV = -M] - G% + Mpfg^'^Szy^k, , Z = Z\ = -3A + ^i? + Mp^g'^S^yuj ■ (13) 

The foliation preserving diffeomorphism includes the three dimensional spatial diffeomorphism as a part of it. As a 
result, Z^j satisfies the generalized Bianchi identity, 

D,Z^^ = 0. (14) 
For convenience, we decompose the spatial metric and the extrinsic curvature as 

= a2(i)e2C(t.5)^^^.(i,f), (15) 

K', = ^K{t,x)S'^+A'^{t,x), (16) 

where we define C{ti^) so that det7 = 1, and a{t) is defined in ea. (pO|) later. The trace part and the traceless part of 
the definition of the extrinsic curvature lead, respectively, to 



and 



7,, = 27,,.^'=, + ^ (^73k^^N'' + i^kdjN'' - ^7.,5fcA^' j . (18) 

The momentum constraint, i.e. the equation of motion for iV*, is 

DjK^^ - XdiK = 0. (19) 
According to the decomposition (1161) . the momentum constraint is rewritten as 



d,A^^ + ZA^^djC - \A\i'''da,k - ^ (3A - 1) d,K = 0. (20) 

As a consistency check, it is instructive to calculate time derivatives of lndet7, 7^^ — 7^^ and ^ikA'^^ — jjkA\ 
without using the fact that they actually vanish. The results are 

d^A\ = ~KA\, 
a±(lndet7) = 2A\, 

^±il^J-JJ^) = 2{-f,kA''^--fjkA\), 

9±(7.fe^', - = -K{l^kA^, ~ ijkA^i) + 2{^uA\ ~ iuA\)A\ - 2(7,7^', - iuA\)A\ 

+2(7« - li,.)A\A'^ + ^(7.fe^^ - lii.A\)d,N' - ^{i,kA\ - 7z.A^)5.iV' 

-\in^kA^, - l,kA\)diN' - iz(7,, - 7,,). (21) 

The right hand side of each equation vanishes when A' ^ , In det 7, 7^ — 7^^ and ^ikA^j — jjuA^ vanish. Therefore, the 
evolution equations we have derived are consistent with vanishing ^'j, In det 7, 7^ — jji and jikA'^, — jj^A^. 
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III. GRADIENT EXPANSION 



The gradient expansion is the method to analyze the full non-linear dynamics at large scale. In the gradient 
expansion, we consider nonlinear perturbation around a flat Friedmann-Robertson- Walker background and suppose 
that the characteristic spatial scale L of the perturbation is much larger than the Hubble horizon scale \/H . To make 
the argument transparent, we introduce a small parameter e defined by e 1/{HL) and expand all relevant quantities 
and equations with respect to e. For example, a spatial derivative acted on a relevant quantity raises the order of e 
and thus is counted as 0{e). We then solve the equations order by order in gradient expansion. 



A. Gauge fixing 

The foliation preserving diffeomorphism invariance is, like all other gauge symmetries, redundancy of descriptions. In 
order to extract physical quantities and statements, we thus need to eliminate gauge freedom by imposing appropriate 
gauge condition. In this paper we adopt the synchronous gauge, or the Gaussian normal coordinate system, by setting 
the lapse to unity and the shift to zero. 

N =1, N' = 0. (22) 

This fixes the time coordinates but does not completely fix the spatial coordinates. There still remains gauge freedom 
of time-independent spatial diffeomorphism, corresponding to the change of coordinates on the initial constant-time 
hypersurface. This residual gauge degree of freedom will be discussed later. 
In this gauge our basic equations (|9]), ([T0|) . (fT7|) and (|18|) are simplified as 

(3A - l)dtK = -i(3A - 1)K' - ^A^^A\ (23) 

dtA'^ = -KA'^+Z'j-^ZS'^, (24) 

dtC = -^ + \K, (25) 
a 6 

^tl^J = 27,feA*^^-. (26) 

Hereafter we assume that A ^ 1/3. Actually, as already explained in the introduction, the regime of physical interest 
is A > 1. 



B. Order analysis 



In order to expand the equations and to write down equations in each order of gradient expansion, we need to know 
the orders of all relevant variables. Therefore, we begin with the order analysis to determine them. 

Since we are interested in the spacetime which is not so much different from the exact Friedmann universe, we 
suppose that 

^a^J = 0(e). (27) 

Substituting this into ea.(|26p. we obtain 

= 0(e). (28) 

Then, the constraint equation (|20p implies that 

d^K = 0(e2). (29) 
In other words, K'^^^ depends on t only. This fact enables us to define a{t) by 

3^^ = X(o)(=3i/(t)). (30) 
a{t) 

With this definition of a{t), eq. ([25)) leads to 

dtC = 0(e). (31) 
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In summary we have the following expansion. 

C = C^"^ {^) + {t, x) + e^C^^) {t, f ) + • • • , (32) 

7., = f.j {x) + ctI;^ (i, x) + e^jlf {t,x) + --- , (33) 

K ^ ■iH{t) + eK^^\t,x)+e^K^^\t,x) + --- , (34) 

A', = eA^^' V (i, X ) + £2^(2) \^{t,x) + --- , (35) 

where a quantity with the upper index (n) is rt-th order in gradient expansion. 



C. Equations in each order 

We have found the orders of all physical quantities. Substituting this into the evolution equations (|23ll26p . we can 
obtain the evolution equations in each order. 

In the zero-th order of gradient expansion we have 

{■i\~\)(dtH+lHA=K. (36) 



The first integral of this equation leads to 



2 



where C is an integration constant. The second term in the right hand side of this equation is the "dark matter as 
an integration constant" [5|. 

The n-th (n > 1) order equations are written as 



V / 2^ 2 3A-1^ ^ ' 3A-1' ^ ' 

n — 1 ^ 



p=i 



(40) 



n-l 

d,^^ = 2 5]7lr^A("-^)V (41) 

p=0 

In a similar way, from eq. (|20p we obtain the {n + l)-th (n > 1) order momentum constraint equation as 

5,a(")^', + 3^a('')^9,C("-^') - A(^)^",(7"i)(«)"=a,7]r''^'^ - i(3A_ l)a,if(") ^ 0, (42) 

p—1 p—1 9—0 

where (7~i)(")*-' is the n-th order part of the inverse of 7^^-, i.e. the inverse (7"^)'-' is expanded as 

(7-1)*^" =/*^"+e(7-i)(i)^J +e2(7-i)(2)»J +... , (43) 

where = (^-i)^*^ is the inverse of fij. It is easy to show that (7 ^)'^"-' {n>\) satisfies the following differential 
equation. 

n 

p=i 

There are some useful identities. The generalized Bianchi identity (|14p leads to 

p=0 ^ ^ p=0 g=0 
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where z'"-**^ and Z^") are the n-th order parts of and respectively. The conditions j4\ = 0, (?ilndet7 = 0, 
"iij — = 0, jikA'^j — "ijkA^i = and — "fjkA''i["f^^)'-'^ — lead to the following identities. 

p=0 

n— 1 n~ln—p—l 

^(^(P)^(n-P)fe_^(p)^(n-P)fe^)^0^ ^(")^_^ ^ 7g)A(™'\(7-i)(^)'^ = 0. (46) 

p=0 p=0 g=0 



D. 0(e) solution 

For the first order (n = f), eas. pSMTj) are reduced to 

dt (a^Xd)) = 0, (47) 

at(a3^(i)v) = 0, (48) 

dtC^'^ = (49) 

a,7« = 2/,,A(i)V (50) 



Note that Z^^-*' = 0. Integrating these equations, we obtain 



„(i) _ Of. /"* , .,(1) 



7^;^ - 2/,,C-- / (54) 



where the integration constants C^^\ C^^^^j, Cln Tinlj depend on the spatial coordinates a;* only, and satisfy 

d''>\ = 0, Ud'^)^f,kd'^',. (55) 

The two integration constants, (^|^^ and Ji^ij^ can be absorbed into the zero-th order counterparts, Cin ^^'^ linij- 
Thus, without loss of generality, we can set 

C£^-0, 7l;ll,=0. (56) 

Finally, the momentum constraint equation (j42p with n — 1 leads to the following relation among the remaining 
integration constants, C^^', C^^^^j, C^"^ ^^^d /y . 

+ 3C(^^^,,9,C^°^ - {c^'^'ifd^ljk - ^ (3A - 1) = 0, (57) 

where is the inverse of fij . 
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E. O(e^) solution 



With the first order solution obtained in the previous subsection, we can solve the second order equations. In the 
second order, eas. pSMTt become 



2^ ' 2((3A-1)) ' 



1 



2(3A- 1) 



(58) 

(59) 

(60) 
(61) 



where ^ and R are Ricci tensor and Ricci scalar constructed from the 0-th order conformally-transformed metric 
a~^glj^ = e^""' ^ fij, and we have used the fact that z''^'' ' ^ ~ — a^^ (^'j ~ integrating these equations 

we obtain 

J (2) ' _ 



R 



dt' 



t 



a{t')dt' \ , 



where we have set 



3A- 1 



t' 



a{t")dt" 



a{t")dt" 



^ ^(2) 



„(2) 



= 



(62) 
(63) 
(64) 

(,65) 
(66) 



by redefinition of C^^\ C j, and fij, respectively. 

Provided that the redefined integration constants (C'^^-', C^^''*^, /y) satisfy ((57| up to 0(e'^), one can show 
that the solution (|62ll65p automatically satisfies the third-order momentum constraint equation 

-iA(^)V''=a,:7« - lA^'^\h-')^'^''da,, ~ \A'^'\f^dd,. - \ (3A- l)a.if(^) = 0. (67) 
The proof is given in Appendix [Xl 



F. 0(e") solution (n > 2) 



For general n (> 2), the solution to eqs. p8ll4T|) is 



(n)i 
J 

(n) 
7 , 



a3(t) 



dt'a^it') 



^ n—1 



n-l 



-P) 



dt'a^t') 



n-l 



2(3A- 1) 



(p)i A{n-p)j 



3A- 1 



3 ^ 



di'is:^'^ 



„t n-l 



(p) 4("~p) 



(68) 
(69) 
(70) 
(71) 
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where we have set 



A 



(n)j 



t=ti„ 



(") 



t=tin 



by redefinition of C'^-', C*'^'' 'j, C*-^'' and /ij, respectively. Note that, in subsection IIII E| we have already set 



c 



(1) 



t=ti„ 



(72) 



(73) 



by redefinition of C'-^'' and fij, respectively. 

The initial condition for 7^"-' (n > 1) implies that 7ij|(^(. — fij, that (7^^)'-'|j_j = and that 
(-y-i)(")«i| 



ij^j = (n > 1). Therefore, for n > 1, the solution to is 

/ dt'^A^P)\(7-i)("-f)^^ 



(74) 



As shown in Appendix 1X1 the solution (|68II71I) automatically satisfies the (n + l)-th order momentum constraint 
equation P^ . provided that the redefined integration constants {C^^\ c'^^* , fij) satisfy ([57| up to 0(e"+^). 



G. Number of physical degrees of freedom 

Our solution involves functions (^'^\x), fij{x), C^^^(x) and C^^''*j(a;) of spatial coordinates as integration 'constant'. 
They are subject to the constraint ((57)) . Also, as stated just after (P^ . our gauge fixing condition ([2^ leaves time- 
independent spatial diffeomorphism as residual gauge freedom. Therefore, the number of physical degrees of freedom 
included in each integration 'constant' is 

1 scalar growing mode — 1 component , 

2 tensor growing modes — 5 components — 3 gauge , 

1 scalar decaying mode = 1 component , 

2 tensor decaying modes = 5 components — 3 constraints . (75) 

This is consistent with the fact that the Hofava-Lifshitz gravity includes not only a tensor graviton (2 propagating 
degrees of freedom) but also a scalar graviton (1 propagating degree of freedom). 



C(°)(f) 

ftjix) 

c(°)(f) 



IV. CONSERVED CURVATURE PERTURBATION 



In this section we propose a definition of nonlinear curvature perturbation and show that it is conserved up to the 
first order in the gradient expansion. This statement holds for any values of A. Because of the conservation, this 
quantity is expected to be useful for the analysis of superhorizon evolution of nonlinear perturbation. 



A. Definition 



In general relativity, it is known that the curvature perturbation in the uniform density slice conserves up to the 
first order in the gradient expansion |7|]. Motivated by this fact, we define the quantity TZ in Hofava-Lifshitz gravity 
by 



7^(^, x) = ({i, x) + In 



a{i) 



ait) 



(76) 



where t{t, x) = t + 0(e) is the solution to 



p'^lit) + Sp,^{i,x)^pZ{t) 



- JO) , 



(77) 
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and 



(78) 



In the following we shall show that TZ is indeed conserved up to the first order in the gradient expansion by explicitly 
calculating it. 

B. Concrete expression and conservation up to 0(e) 
According to the gradient expansion, we expand Spdm, t and TZ as 



ri=0 



where a quantity with the superscript (n) is of order 0(e"). 
At the order 0(e°), we obtain 



and this is constant in time. 
At the order 0{e), we obtain 



.(1) 



Thus, the time derivative of TiP-^ is shown to vanish as 



7^(") = c^")(x). 



(79) 



(80) 



(82) 



where we have used the time derivative of p6p to show the last equality. 

We therefore conclude that TZ defined in (|76p is conserved up to the first order in the gradient expansion. This 
statement holds for any values of A. 



V. SUMMARY AND DISCUSSION 



We have performed a fully nonlinear analysis of superhorizon perturbation in Hofava-Lifshitz gravity by using the 
gradient expansion technique. In § IIIII we have presented a concrete expression for the solution of gravity equations 
up to the second order in the gradient expansion. We have also proven that the solution can be extended to any order 
of the gradient expansion, by showing that the solution to the dynamical equation satisfies the constraint equation at 
each order. 

Based on the result, in §[lV]we have proposed a definition of nonlinear curvature perturbation TZ in Hofava-Lifshitz 
gravity and have shown that it is conserved up to the first order in the gradient expansion. 

It is known that, in the limit A — > 1, the scalar graviton gets strongly coupled and the usual metric perturbation 
breaks down in the scalar graviton sector Q. Here, we stress that the breakdown of the perturbative expansion does 
not necessarily lead to a loss of predictability since all coefficients of infinite number of terms in the perturbative 
expansion can be written in terms of finite parameters in the action if the theory is renormalizable. Indeed, for 
spherically symmetric, static, vacuum configurations, one can perform fully nonlinear analysis to show that the limit 
A — ^ 1 is continuous and the general relativity is recovered in the limit |^]. This result may be considered as an 
analogue of Vainshtein effect and suggests the possibility that the scalar graviton may safely be decoupled from the 
rest of the world, i.e. the tensor graviton and the matter sector, in the limit. 

The result of the present paper is based on the fully nonlinear analysis and may be considered as yet another 
example of the analogue of Vainshtein effect. Up to any order of the gradient expansion, the equations of motion and 
their solutions are manifestly regular in the limit A ^ 1. The solutions reduce to those in general relativity coupled 
to dark matter in the limit at low energy. 
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In the present paper we have concentrated on the pure gravity system in the projectable Hofava-Lifshitz theory. 
Because of the existence of the scalar graviton, this simple system is still rich enough as a testing ground for the 
analogue of Vainshtein effect. Indeed, the so called "dark matter as integration constant" Q drives non-trivial 
cosmological dynamics in this system, and thus the nonlinear analysis presented in the present paper provides a 
convincing evidence for the analogue of Vainshtein effect. It is certainly interesting and important to extend the 
nonlinear analysis to more general situations with matter contents [lol |. 
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Appendix A: (n + l)-th order momentum constraint (n > 2) 

In this appendix, by induction we prove that the n-th order solution (j68ll7ip satisfies the (nH-l)-th order momentum 
constraint equation (|42l) for n > 2. 

The basic logic of the proof is to rewrite the left hand side of the {n + l)-th order constraint (H^ as a linear 
combination of lower order constraints by using the explicit solution (|68ll7ip . For this purpose we shall use the 
generalized Bianchi identity (|^5|) and other identities We also use the following identity for functions f{t) and 

g{t) satisfying a^{Un)f{tin)9{Un) = 0: 



/WffW = ^ / dt'a^it') [a{tr'dAa\t')f{t')) ■ g{t') + f{t') ■ dt^git')] . (Al) 

By applying the identity m to (,/(*),. gW) = a.C^"-^')) and {f{t),g{t)) = (A^") ^, (7-i)(*) "=9,7]^"'""'), 

the left hand side of the {n + l)-th order momcntmn constraint equation ()42|) is rewritten as 

n n—p 



" -| It ^ ^ 

p—1 p—1 q—0 

^ ft ( n 

-3^^ dt'aHt') 1 3^^ [a~^dt' (a^A^^^^) d,&-P^ + A^^^'^d, [dt>&-P^ 

n n—p 



(A2) 
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Using (|68ll69p . (|40ll4ip and (|44)) . this is further rewritten as 



c. 



(n+1) 



n / p-1 \ n-1 y V 

.p=2 \ 9=1 / p=l ^ ^ 



n n— p 

^EE 

p=l q=0 



p-1 



(n-p-q) 



7i—p—q—l 



\ r=l 

-i(3A- i)a, pg^'^^^^"^"^) - (-E 

n n— p ~j 



p=l q=0 J 

where we have used the generahzed Bianchi identity (|45p . By using the identities (|46p we finally obtain 



C. 



(n+l) 



1 /■'^ 



Since we already know that C,^^'' = under the condition ((57|) . this is enough to prove C^"+^^ = for n > 2. 



(A3) 



(A4) 
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